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We analyze the optical spring characteristics of a double pumped Fabry-Perot cavity. A double- 
resonance optical spring occurs when the optical spring frequency and the detuning frequency of 
the cavity coincide. We formulate a simple criterion for the stability of an optical spring and apply 
it to the double resonance regime. Double resonance configurations are very promising for future 
gravitational wave detectors as they allow us to surpass the Standard Quantum Limit. We show 
that stable double resonance can be demonstrated in middle scale prototype interferometers such 
as the Glasgow 10m- Prototype, Gingin High Optical Power Test Facility or the AEI 10m Prototype 
Interferometer before being implemented in future gravitational wave detectors. 
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I. INTRODUCTION 

Currently the search for gravitational radiation from 
astrophysical sources is conducted with the first- 
generation Earth-based laser interferometers (LIGO in 
USA [1-3], VIRGO in Italy [4, 5], GEO-600 in Germany 
[6, 7], TAMA-300 in Japan [8, 9] and ACIGA in Australia 
[10, 11]). The development of the second-generation GW 
detectors (Advanced LIGO [12, 13], Advanced Virgo [14], 
GEO-HF [15] and LCGT [16]) is well underway. 

The sensitivity of the first-generation detectors is lim- 
ited by noises sources of various nature: seismic and sus- 
pension thermal noise noise at low frequencies (below 
~ 50 Hz), thermal noise in suspensions, bulks and coat- 
ings of the mirrors (~ 50 — 200 Hz), photon shot noise 
(above ~ 200 Hz). It is expected that the sensitivity of 
the second-generation detectors will be ultimately limited 
by the noise of quantum nature arising due to Heisen- 
berg's uncertainty principle [17-20] over most of the fre- 
quency range of interest. The optimum between mea- 
surement noise (photon shot noise) and back-action noise 
(radiation pressure noise) is called the Standard Quan- 
tum Limit (SQL). This level is expected to be reached 
in the forthcoming second generation of large-scale laser- 
interferometric gravitational-wave detectors. Third gen- 
eration detectors, such as the Einstein Telescope [21] aim 
to significantly surpass the SQL over a wide frequency 
range [22]. 

The most promising methods to overcome the SQL rely 
on the implementation of optical (pondcromotivc) rigid- 
ity [20, 23-25] which effectively turns the test masses of 
a gravitational-wave detector into harmonic oscillators 
producing a gain in sensitivity [26-31]. In order to un- 
derstand it note that the formula for the sensitivity £(Q) 
of Advanced LIGO interferometer (see Fig. 1) consists of 
two terms: 
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Here Sh is the total single-sided spectral density of the 



noise recalculated to the strain density h of the gravita- 
tional wave /igQ L = 2H/ \iL 2 n 2 — the value of Sh corre- 
sponds to the SQL sensitivity for the case of free masses, 
fl is the observation frequency, L the arm length of the 
interferometer, fx = m/4 is the reduced mass, m is the 
mass of each mirror of the arm cavities, S x and S F are the 
single-sided spectral densities of the measurement noise 
and the back action noise, respectively and /LtA'(f2) is the 
optical rigidity. The relation between the Fourier trans- 
forms of the position x and the force F can be described 
as - fi 2 + K(VL))x(VL) = F(fi). For simplicity we 
assume here that the measurement and the back action 
noise are not correlated. Then the uncertainty relation 
is given by 
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Substituting (2) in form S F (fl) > h 2 / S x (fl) into (1) we 
can find the minimal value of the sensitivity (after 
optimization over S x ) which is better than the SQL 
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in the bandwidth Afl where \K(Q) — Q 2 \ < f2 2 , i.e. close 
to optical spring resonance. In other words a harmonic 
oscillator provides a gain in sensitivity for near-resonance 
signals, equal to 
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Usually the two resonances are at two separate frequen- 
cies. The gain in sensitivity by means of the optical rigid- 
ity was examined in [28, 29, 32] for the wide-band regime 
under conditions Ail ~ f2 and £ > 1/2. 

The case of so-called double resonance [27, 30, 31], 
when an sophisticated frequency dependence of the op- 
tical rigidity allows us to obtain two close or coinciding 
resonance frequencies and to get better gain in sensitiv- 
ity, is described by the following formula: 
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It is important to note that the oscillatory sensitivity 
gain, described in Equation (4) does not provide any gain 
in the signal to noise ratio in case the signal has a band- 
width larger than Aft (because the signal to noise ratio 
scales as the sensitivity gain multiplied by the bandwidth 
AO). However, for a double resonance configuration with 
a sensitivity described by Equation (5) the signal to noise 
ratio increases with decreasing bandwidth proportional 
to Ail" 1 [27, 30, 31]. Please note that in the Advanced 
LIGO configuration optical rigidity can be created eas- 
ily through microscopic position changes of the signal 
recycling mirror [28-31] as it is already being used in 
GEO 600 detector [6]. 

A single optical spring always causes instability. This 
can qualitatively be explained [23-25] by taking into ac- 
count that the optical rigidity is not introduced instan- 
taneously, but with a delay of the relaxation time r* of 
the optical resonator. Therefore, the evolution of the free 
mass position x with the optical rigidity K may approx- 
imately be described by the equation 



mx{t)+Kx{t-T*) = 0. 



(6) 



Expanding the second term in a series x(t — r*) ~ x(t) — 
r*x(t) we can rewrite the previous equation in the form: 



mx(t)- Kt* x(t) + K x(t) = 0. 



(7) 



Obviously, the term —Kt* x corresponds to a negative 
damping force. Please note that in the opposite case 
of negative optical rigidity a positive damping force is 
introduced (the sign of the rigidity correlates with the 
sign of the detuning). 

The instability (negative damping) can be compen- 
sated by incorporating a linear feedback control loop and 
in the ideal case (no additional noise is introduced by 
the feedback) it would not modify the noise spectrum of 
a GW detector [29]. In practice, however, the need for 
control gain at frequencies inside the detection band can 
cause undesirable complexity in the control system or can 
introduce additional classical noise. 

An alternative way to suppress the instability was pro- 
posed [33] and experimentally demonstrated [34], by in- 
jecting a second carrier field from the bright port (see 
Fig. 1) in order to create a relatively small additional neg- 
ative rigidity component, thus the total rigidity (of both 
lasers together) remains positive, but at the same time 
introduce a relatively large additional positive damping 
component to make total damping positive. The main 
purpose of the second carrier is to create a second op- 
tical spring that forms a stable optical spring together 
with the first one — even though each individual optical 
spring, acting alone, would be unstable. Both carriers are 
assumed to have different polarizations, so that there is 
no direct coupling between the two fields (although they 
both directly couple to the mirrors). 

In this paper we further analyze the stability of an op- 
tical spring created by a double pump. The regime of 
stable double resonance analyzed in this paper provides 




FIG. 1: Scheme of an Advanced LIGO interferometer pumped 
by two lasers. The main laser is detuned to give a positive 
optical rigidity and negative damping (it is tuned on the right 
slope of resonance curve) while the auxiliary laser is detuned 
to give negative optical rigidity and positive damping. 



the possibility to decrease the bandwidth of the resonance 
curve dramatically and, hence, to increase the signal to 
noise ratio even for signals with a wide bandwidth (as 
formula (5) predicts). In Sec. II we formulate the simple 
criterion of stability and we apply it to the regimes of sta- 
ble double resonance in Sec. III. In addition we discuss 
the possibilities to observe these regimes using the ex- 
perimental setups of the Gingin and Glasgow Prototype 
interferometers. 



II. STABLE OPTICAL SPRING 

In the following we use notations similar to the ones 
introduced in [28, 35]. We consider an Advanced LIGO 
interferometer with a signal recycling mirror (SRM) hav- 
ing an amplitude reflectivity p and a power recycling 
mirror (PRM) as shown in Fig.l. We assume the mir- 
rors to be without any optical losses. In addition, we 
suppose that both Fabry-Perot (FP) cavities in the east 
and north arms are identical, each input mirror featur- 
ing an amplitude transmittance T <C 1 and an reflectiv- 
ity R = Vl — T 2 , while each end mirror is completely 
reflecting. The end and input mirrors have an identi- 
cal mass m and in absence of the laser pumps they can 
move as free masses. The PRM is used only to increase 
the average power incident onto the beam splitter and we 
assume that no fluctuational fields from the laser (west 
arm) reach the detector in the south arm. The mean 
frequency of each pump Wx,2 is equal to one of the eigen- 
frequencies of the FP cavities. L is the distance between 
the mirrors in the arms (4 km for (Advanced) LIGO), I is 
the mean distance between the SR mirror and the beam 
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splitter. Wc also introduce the following notations: 



m(T\ + A 2 )/8 is a convenient multilpier) : 
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Here 70 is the relaxation rate of a single FP cavity, 
A12 are the dctunings introduced by the SR mirror, 
Tx2 are the relaxation rates of the differential modes of 
the interferometer for each pump. The frequency of the 
differential mode depends on the arm length difference 
z = xe — He — ( x n — Un) an d its bandwidth and de- 
tuning are controlled by position of SR mirror, which is 
described by phase advance 01,2 arising from differential 
mode detuning (0i,2 = corresponds to resonance). We 
assume that 0i.2 do not depend on f2 due to the small 
length I of the SR cavity: I -C L. It is worthwhile not- 
ing that the detunings Ai,2 and relaxation rates Ti,2 can 
differ for different pumps. 

The mechanical evolution of the differential arm length 
z = xe — He — xn + Un (see notations on Fig. 1) in the 
frequency domain is described by 

z=^(fl)(F 1 + F 2 +F s ), F s =fin 2 Lh, (12a) 
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where h is a dimensionless metric of a gravitational wave, 
\& is the mechanical susceptibility (-Ki,2 are the frequency 
dependent spring coefficients created by the correspond- 
ing laser pump); ix,2 are the average optical arm cavity 
powers of the carrier 1 or 2 with the mean optical fre- 
quencies wi, 2 ; Fx, F 2 are the back action forces caused 
by the first and second pump. 

The eigen frequencies are solutions of the characteristic 
equation 
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Roots corresponding to stable oscillations must have a 
negative imaginary part (as we assume the time depen- 
dence of position to be ~ e~ lUt ). 

The equation (13) can be rewritten using the di- 
mensionless susceptibility \ — ^^(r 2 + A 2 )/8 (here 
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(14b) 
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We choose the dimensionless frequency x so that for a 
single pump (i.e. Q = 0) the double (unstable) resonance 
would take place at x 2 ~ 1 if P = 1 [30] (in formal limit 
j^O when susceptibility is a pure real value). 




FIG. 2: Roots correspond to a stable (slightly damped) oscil- 
lation if and Yj have the same signs near roots of Y r , i.e. 
plots of Y r and Yi correspond to each other as shown on plot. 

In order to find the eigen-frequencies one has to solve 
the equation Y(x) = and to analyze for which condi- 
tions the roots are stable. Undoubtedly the roots of this 
equation can easily be obtained by applying numerical 
methods. However, starting from a set of input parame- 
ters (powers and detunings of each pump) the numerical 
solution will provide a set of frequencies with arbitrary 
imaginary parts. The problem is to identify the values of 
the system parameters that actually provide stable roots 
with negative imaginary parts. In the following we for- 
mulate a criterion which allows us in a quite easy way 
to properly estimate a set of parameters fulfilling our re- 
quirements. 

In the general case, Y r is a third order polynomial of 
x 2 , so all of it's roots (which we assume to be real) can 
be written in the form of (j = 1, 2, 3). An example 

of Y r as a function of x 2 is plotted on Fig. 2. 

In order to apply a successive approximation method 
we assume that the values of the function Y; arc small 
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in the regions close to the roots of Y r . If this is the 
case, then the roots of the equation Y(x) = can be 
found using the successive approximations method taking 
6 roots ±2y of the equation Y r (x) = as a zeroth order 
approximation. In the first order approximation we add 
a small imaginary part to each root: 



xy + iSj , xf 1 



(0) 



(15) 



After substituting this into Eq. (14b) we find 5ji 



(0) 



= 0, 



Y (x {oy 



(16) 



where Y r ' (x 



( o)\ = dY r 
j I ~ d{x 2 ) 



< x fT 



In order to fulfill the condition for a stable oscillation, Sj 
have to be negative as we required in (16) (recall, as we 
assume the time dependence of position to be ~ e~* n *, 
hence, eigen frequency (xj) corresponds to oscillations 
~ e~ lXjt and negative imaginary part means damping). 
Therefore, stable roots with Sj < will be achieved in 
the case when the functions Y r ' and Y have the same 
signs in the regions close to roots of Y r , i.e. their plots 
have to relate to each other as shown on Fig. 2. 

We would like to emphasize that the obtained crite- 
rion is very simple and that it allows us to estimate 
whether a certain set of input parameters can provide 
roots (complex frequencies) with imaginary parts of the 
desired signs without the need to solve the character- 
istic equation. The proposed criterion is very close to 
the Raus-Gurwitz criteria [36], however, it is more obvi- 
ous and more useful for the particular analysis of optical 
spring stability. 

The stability criterion described in Eq. (16) can easily 
be generalized for a larger number of optical pumps. For 
example, for three pumps the functions Y r and Y are 
parabolas of fourth and third order correspondingly (in 
respect to x 2 ) and the roots will be stable if the functions 
Y r ' and Y have the same sign at the regions close to the 
roots of Y r as shown in Fig. 3. 




FIG. 3: For the case of 3 pumps the roots are stable if in ac- 
cordance with criterion (16) the plots of Y r and Yi correspond 
to each other as shown. 

Obviously, the regimes of stable optical spring consid- 
ered in [33] fulfil the condition described in Eq. (16). 





FIG. 4: "Ideal" double resonance: curve Y r toughes z-axis in 
point zq: a) left root z\ of Y coincides with zq; b) curve Y 
touches z-axis in the same point zq. 



However, our criterion might be useful to identify config- 
urations that can be experimentally implemented avoid- 
ing any feedback or for other interesting regimes such as 
double resonance or negative inertia [37]. 

TABLE I: Parameters of Advanced LIGO used in this paper. 



Parameter 


Value 


Arm length, L 


4 km 


Mass of each mirror, m 


40 kg 


ITM amplitude transmittance, T 


V0.005 


Relaxation rate of single FP cavity, 70 (8) 


94 s _1 


SRM amplitude reflectivity p 


v / 093 


Optical wavelength, A 


1064 nm 



III. STABLE DOUBLE RESONANCE 

Oscillators in the double resonance regime have inter- 
esting properties. The resonant force, with the time de- 
pendence f = Fq cos f2 rcs £, acting on such an oscillator 
produces a displacement z^i = — Foi 2 /(8m) cos tt lcs t ~ 
i 2 , which is much larger than the resonant displacement 
2 co = — F t/(2m) sinf2 rcs ~ t of a conventional oscilla- 
tor under action of the same force [30]. In the frequency 
domain the resonant gain of a conventional oscillator is 
~ f2 res /Af2 (i.e. inverse proportional to the bandwidth 
AD, of the resonance curve) whereas the resonant gain 
of an oscillator in the double resonance regime is much 

larger: ~ (fi res /AQ) (i.e. inverse proportional to the 

square of the bandwidth (Afi) 2 ). This feature, as men- 
tioned in the introduction, is responsible for the wide 
bandwidth gain in signal to noise ratio. 

The double resonance takes place if the susceptibil- 
ity has a pole of second order [27, 30]. In terms of our 
stability criterion this means for the case of the double 
resonance that the plot of Y r as function of x 2 touches 
the horizontal axis at the resonance frequency x$ and 
Y is equal to zero at this point. Fig. 4 illustrates two 
possible ways of behaviour of Y: in Fig. 4a and 4b Y 
crosses and touches the horizontal axis at x 2 = Xq, re- 
spectively. However, it is difficult to exactly realize the 
regimes shown in Fig. 4 — some discrepancies like shown 
on Fig. 2 are inevitable. 
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In general, double resonance takes place when any two 
of three roots [x^] 2 , [a^] 2 , [ x 3°^] 2 coincide. In this 
section we consider the particular case of close smaller 
roots , [z^] °f Y r (our preliminary analysis 

shows that the other case of close roots [x 2 ' ) ] 2 , f 2 -^] 2 
can not provide a stable set of roots): 
2 



„(0) 
'1,2 



xi±e, e < 1 



(17) 



where Xq is a middle point between the roots 
Hence, Y r may be presented in the form 



(o) 

C l,2 



{x 2 — Xq - e)(x 2 — x 2 + e) 



(18) 




FIG. 5: Detailed plot of two close roots of Y r (14c). The roots 
of Yi (14d) have to be shifted to the right relative to the roots 
of Y r in order to fulfill the stability criterion (16). 

In order to fulfill the criterion (16) one has to arrange 
the parabolas along the x 2 axis in such a way that the 
roots of Yi would be placed close to roots of Y r (see 
Fig. 5). Shifting Yi along the vertical axis allows us to 
manipulate the damping at the resonant frequencies, be- 
cause the values of Yi at the roots of Y r dehne the me- 
chanical damping of the oscillator (16). The problem of 
properly arranging this two curves is solvable as there are 
enough parameters for manipulation available: the two 
detunings and the two powers of the two pumps (i.e. the 
dimensionless parameters P, Q, r/, v). An example of a 
step by step calculation is presented in Appendix A. For 
the following dimensionless parameters 

P ~ 0.77, Q ~ -0.027, g = 0.6, d~ 1.39, (19a) 

?7~0.24, i/ = 0.5, e = 0.05. (19b) 

we get the eigen frequencies: 



n 3 



532.3 s" 
482.1s" 
961.8 s" 



Si 
62 



16.8 s" J 
15.4 s" 1 
314.2 s" 



(20) 
(21) 
(22) 



Using Advanced LIGO parameters as listed in Table I 
the set (19) can be recalculated, giving the following re- 
sults: 



Ti ~ 230s" 1 , 
T 2 ~ 115s" 1 , 
(0i 1.49, 



Ai 
A 2 



1063s" 
-526 s" 
-1.39), 



h : 
h 



863 kW, 
i 61 kW, 




FIG. 6: Trace 1: susceptibility |^(/)| as function of frequency 
/ = Q./2-k for the parameters given in Eqn (19). Traces 2,3: 
same as trace 1 but with e — > e/2 and with e — > e/4. Trace 
4: susceptibility of free mass (4/m(27r/) 2 ). 



Fig. 6 shows the corresponding susceptibility vfi" 1 as a 
function of frequency / = CI/2t: for the listed parame- 
ters. It is worth emphasising that it is possible to realize a 
stable double resonance with arbitrary small damping by 
just choosing smaller values for e. Fig. 6 includes the sus- 
ceptibility for several values of parameter e: e— > e/2, e/4, 
as well as the susceptibility of the free masses as compar- 
ison. 



A. Possibilities of experimental observation 

It would be interesting to demonstrate and observe a 
stable double resonance in an experiment. For example, 
this could be done in Gingin High Optical Power Test 
Facility [38]. The initial formulas (14) are valid for the 
Gingin topology with one exception. Only one arm of the 
Michclson interferometer (with Fabry- Perot arm cavities) 
is used in Gingin, hence, one has to define the power 
parameters P and Q a factor 2 smaller: 



16fc]I]Ai 



mL(T 2 



A?) 



Q = 



16fc 2 / 2 A 



mL(T\ 



A?) 



(23) 



The Gingin facility features test masses of m = 0.8 kg and 
arm length is L = 80 m. Then, assuming Ti = 400 s _1 , 
it is easy to recalculate the parameters (19) for Gingin: 



Ti 

r 2 



400 s" 
200 s" 



Ai 
A 2 



843 s" 1 , h 
-911s" 1 , Jo 



3.59 kW, (24) 
: 0.254 kW. (25) 



It seems that these parameters may be relatively easy 
realized using the Gingin facility. 
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set of experimental parameters: 



FIG. 7: Scheme of the optical rigidity experiment within the 
Glasgow Prototype Interferometer. Only the end mirror can 
move, while the mirrors PRM and IM can be considered as 
unmovable, because their masses are 30 times larger than that 
of the end mirror. 



TABLE II: Parameters of Glasgow Prototype Interferometer. 



Parameter 


Value 


Arm length, L 


10 m 


Distance between PRM and ITM, I 


5 m 


Mass of movable mirror, m 


100 g 


Masses of PRM and IM, M 


3 kg 


PRM amplitude transmittance Tprm 


VoM 


ITM amplitude transmittance, Tim 




Input power, /input 


> i w, 


Optical wavelength, A 


1064 nm 



Another experimental setup that can be used to ob- 
serve a stable double resonance is the Glasgow 10m Pro- 
totype [39, 40]. In contrast to the balanced scheme of 
a full Michclson interferometer such as Advanced LIGO, 
where it is possible to pump the interferometer with the 
symmetric mode tuned to resonance, while observing the 
optical rigidity in the anti- symmetric mode, in the Glas- 
gow Prototype (shown on Fig. 7) the same mode will 
have to be used to pump and to observe rigidity. Re- 
call that the smaller the detuning the larger is the cir- 
culating power, whereas optical rigidity requires a large 
detuning. However, it is possible to set the first cavity 
between PRM and IM to anti-resonance (the relaxation 
rates ri = T 2 = Tpr,m7imc/4L reach their minimum) 
and to create the detuning by shifting the position of 
the end mirror. In this situation the equations (14) can 
be used with v = 1. Taking into account that practically 
only the light end mirror can move the power parameters 
P, Q have to be redefined, reduced by factor 4: 



P = 



8fei/iAi 



mL{T\ 



A?) 



Q = 



8fc 2 / 2 A 2 



mL{T\ 



(26a) 



Then we can apply the scheme presented in Appendix A. 
In particular the stable double resonance is characterised 
by the following set of dimensionless parameters: 



P ~ 2.43, Q ~ -0.51, g = 0.625, d 
77 ~ 0.71, v=l, e = 0.01. 



1.27, 



(27a) 
(27b) 



Using the actual parameters of the Glasgow Prototype, 
as presented in Table II, we can calculate the following 



ri = r 2 - 

■^1 input 



937 s"\ 
- 1.28 W, 



A 2 ~ 4138s" 1 , A 2 



'2 input 



0.33 W. 



-3450s" 1 , 
(28) 



are the input powers of the two 



Here I\ input) I2 input 

pump lasers. 

Unfortunately, the observation of a stable double res- 
onance is more difficult for smaller scale experiments. 
For example, one may consider an experimental setup 
as described in [34] of a short Fabry-Pcrot cavity (with- 
out power recycling) and a light-weight movable mir- 
ror: the power transmittance of the input mirror is 
T 2 = 0.8 x 10~ 3 , the distance between the mirrors is 
L ~ 0.9 m and the mass of the movable mirror is m = 1 g. 
Then the dimensionless parameters (27) can be recalcu- 
lated as follows: 



ri = T 2 ~ 66.7 x 10 d s~ 
Ai ~ 294 x 10 3 s~\ A 2 



-245 x 10V 



h 



input 



2.64 kW, J a ^ ~ 1.26 kW. 



(29) 
(30) 
(31) 



Obviously these input power values are impermissibly 
huge to be realised in an experiment. In order to find 
some regularities concerning the experiment scale we can 
rewrite the formula (26) in the following form using the 
value of input power: 



P = 



CfT 2 /; 



1 input 



L 2 h 



input 



m C 3(T eff /4)V (l + A 2 /r 2 )3 mT e 4 ff 



(32) 

where we substitute T± = T^c/AL. For the Glasgow 
prototype the effective transmittance is equal to T e g = 
7prm7im/4. 

In order to find out whether an experiment is realiz- 
able, we emphasized in expression (32) only the propor- 
tionality to the terms, which significantly change with 
changes of the experiment scale. Decreasing the scale 
of the experiment is usually represented by decreasing 
the mass m of the movable mirror and the distance L 
between the mirrors, which approximately compensate 
each others contribution into P. That means that in or- 
der to achieve the needed value of P and thus to observe 
the desired stable double resonance one needs to keep the 
effective transmittance roughly the same as in a bigger 
scale experiment due to the strong dependence of P on 
T ff. This leads to the need for the circulating power to be 
the same as for the bigger scale experiment, which would 
be pretty difficult due to problems caused by absorption 
induced heating of the small mirror. 



IV. CONCLUSION AND OUTLOOK 

We developed a new criterion for the stability of an 
optical spring and applied it to the double resonance 
scheme. Before such a concept can be realised in a large 
scale gravitational wave detector it would be beneficial 
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to demonstrate it on a prototype scale experiment, such 
as the Glasgow Prototype [39, 40], Gingin High Opti- 
cal Power Test Facility [38] or the AEI 10m Prototype 
Interferometer [41]. 

It is worth noting that our criterion for optical spring 
stability may also be applied to the negative inertia 
regime recently proposed by F.Ya. Khalili and colleagues 
[37] (the term "negative inertia" was proposed earlier 
by H. Mullcr-Ebhardt [42]). This directly follows from 
the fact that the negative inertia scheme formally corre- 
sponds to a double resonance configuration with a reso- 
nance frequency of zero. 

The formulated stability criterion is also valid for grav- 
itational wave detectors pumped by more than two lasers. 
This may allow us to realize stable triple ( or quadro ) reso- 
nance, i.e. regimes when three (or four) eigen-frequencies 
coincide or are close to each other. In this case the re- 
sponse to force at the resonance frequency is greater than 
for a double resonance configuration. In the frequency 
domain the resonant gain for a triple resonance has to 
be proportional to ~ (f2 res /AQJ (i.e. inverse propor- 
tional to the third power of the bandwidth (Ml) ). It 
will allow us to further increase the gain in the signal to 
noise ratio for wide bandwidth signals, mentioned in the 
Introduction, by a factor f2 res / Af2 compared to the dou- 
ble resonance. In this paper no sensitivity analysis for 
gravitational detectors was performed, but we plan to do 
this in a future publication. 

An stable double resonance featuring arbitrarily large 
susceptibility may allow us to use it for the production 
of squeezed light by means of the ponderomotive non- 
linearity [34, 43]. However, the obtained squeezing will 
have specific features [33]: One would have to measure 
two optimal quadratures, corresponding to each carrier, 
and squeezing can only be observed by finding their cor- 
rect combination, because they are correlated (entangled) 
through the mechanical degree of freedom. 

Due to the combination of a uniquely large suscepti- 
bility and a low noise level the stable double (or triple) 
resonance may also be applied to macroscopic quantum 
mechanics experiments with mirrors of relatively small 
mass. In particular, it can be used for the observation 
of quantum entanglement between an oscillator and an 
cm. field [44, 45] or for preparation of the mechanical 
oscillator in a non-gaussian quantum state [46] . 

In addition stable double resonance may be a useful 
instrument in other precision measurements, such as the 
measurements of thermal noise of mirror coatings. In 
case the mass m of one cavity mirror is small enough 
compared to the other one, the equation for the variation 
z of the optical distance inside the cavity may be written 
in the frequency domain as follows: 

m ( - ft 2 + K(Q)) z = F ba - mfl 2 zth , (33) 

where z t h are the thermal fluctuations of the mirror sur- 
face in respect to the mirror's center of mass and is 
the back action force. Homodyne detection of the output 



wave yields the information on z. One can obtain reso- 
nant gain of the signal (in this case it is z t h) inside the 
resonance bandwidth (i.e. for il = O res ± Ml/ 2) and in 
case of a stable double resonance this gain is unusually 
large. 



Acknowledgments 

The authors are grateful for fruitful discussion with 
Farid Khalili, Stefan Danilishin and Kenneth Strain. 
A.R's and S.V.'s research has been supported by the 
Russian Foundation for Basic Research Grant No. 08- 
02-00580-a and NSF grant PHY-0651036. S.H. is sup- 
ported by the Science and Technology Facilities Council 
(STFC). 



Appendix A: Details of the double resonance 
stability analysis 

In this appendix we present the example estimation 
of a set of input parameter values needed to realize the 
regime of a stable double resonance. We will deal mainly 
with the functions Y r and Yi in this section. As it was 
mentioned in Sec. Ill we need to arrange these parabolas 
along the rc-axis to achieve close stable frequencies and 
then arrange the parabola Yi along the y-axis to obtain 
the desired values of the imaginary parts of the frequen- 
cies or, other words, the desired bandwidth of the double 
resonance. 

First we assume the two smaller roots of Y r to be close 
enough to each other. That means that the function Y r 
can be expressed in form (18) 

Y r = (z - z Q - e) (z - Zo + e) (z - z 3 ), z = x 2 , (Al) 

where 2e is the distance between the close roots (and e 
is small). We also introduce the parameter d — z 3 — z 
which is defined as the distance between the center of the 
two close roots and the "lonely" root. 

Proper substitution transforms this cubic function into 
the reduced form Y r = y 3 + py + q. In our case the 
substitution is z = y + (2zo + z^)/3. The particular case 
of this equation with z = zo and y = yo (yo is the middle 
of the close roots in terms of the shifted variable) leads 
us to the following equation: 



After the substitution is done the parameters p and q 
are expressed in terms of e and d as follows: 

P=-\d?-e 2 , g = -| d (d 2 -9e 2 ). (A3) 

On the other hand the function Y r is expressed (14c) in 
terms of the input parameters such as the pump powers 
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and the detunings (parameters P, Q, g, rj 1 v). Substi- 
tution of 



y = z + 



2(1 + ri) + isg 2 



transforms the function into the familiar form Y r = 
y + py + q. However, this time the coefheients p and 
q are expressed in terms of the above mentioned input 
parameters. 

Equating the pairs of expressions for p and q (expressed 
in different terms) one can obtain the solutions for the 
pump powers P and Q as functions of the remaining in- 
put parameters as well as e and d. 

At this step the ^-arrangement of the function Y r is 
completed. Indeed, if we set a certain value of d we can 
expect the middle of the close roots (in the shifted vari- 
ables) to be equal to — d/3. Setting a certain value of e 
results in the close roots to be separated by 2e. 




FIG. 8: Detailed plot of the two close roots of Y r (14c): 21,2 = 
zq ± e. The roots of Yi (14d) are separated by the same 
distance 2e, but the trace of Yi itself is shifted by the distance 
ti- 



lt is obvious that in order to obtain a single resonance 
peak, one should choose a distance between the close fre- 
quencies smaller than the mean imaginary part of these 
roots Sj. Since z = x 2 , the distance between the close 
frequencies can be estimated as x 2 — x\ ~ e/ \/zq. Equan- 
tion (16) gives an estimate for 5. Summing this up one 
can write 









5^1,2) 




2*7(zi )2 ) 




4ed 






4e 


2 d 




5i y 


/So ' 



(A5) 



For us regimes with infinitely close frequencies, i.e. 
with e — ¥ 0, are of interest. Along with decreasing the 
distance between frequencies we can also decrease the 
bandwidth of the double resonance or in other words the 
values of Sj, which is possible by decreasing the values 
of Yi. Equation (A5) is the only condition these values 
should obey. 

To ensure this we apply the following conditions to the 
values of Yr. 



Y l (z 1 ) = si 



4e 2 d 



Yi(z 2 ) = -S2 



Ae 2 d 



(A6) 



where the variables s\ and S2 are subject to fitting and 
should be approximately equal to unity. 

This is the final step, resulting in a set of three stable 
frequencies, two of which form the resonant peak of the 
desired bandwidth in the susceptibility. 

It is worth noting that the parameter e is proportional 
to distance the between the roots of Y r and, hence, e 
is proportional to the bandwidth of the resonance. The 
parameters g, v are free and can be used to vary the res- 
onance frequencies or the values of the pumping powers. 



The next step is to arrange Yi thus it provides stability 
of the complex frequencies. According to the stability 
criterion (16) frequencies are expected to be stable if the 
roots of Y r (z\, z 2 ) and Yi {z\, z%) are located in the 
same way as is shown on Fig. 8: 



Z\ < Zi, z 2 < z 2 - 



(A4) 



The values z\ >2 can be estimated using Cardano's formu- 
lae (see Appendix B) and 21,2 are the roots of a quadratic 
equation: 



21,2 



g^V3 t 2(1 + 77) + i/.9 2 
2. 



V + T) 

21,2 = -j— T 

1 + V 



v + rj 
l + i/ 



vP + Q 



After properly arranging Yi thus it complies with the con- 
dition (A4) the second step is completed. We obtained 
three stable frequencies, two of which are separated by 
the desired distance of 2e. 



Appendix B: Cubic equation: Cardano's formula 



The cubic equation 



x A + ax + bx + c = 



(Bl) 



(B2) 



may be written in the depressed cubic form by substitu- 
tion of x — y — a/3 [36]: 

= y 3 + py + q, 

a 2 /a\ 3 ab , 

P=-y+6, * = 2(-) - y+c (B3) 

The roots of this depressed cubic equation are given by: 
A+B r A-B 

-^—±iV3^—, y 3 =A + B, (B4a) 



2/1,2 



where 



1/3 



<1 



\-^d) , (B4b) 



D 



(f) 



3 , qs 2 



(B4c) 
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For A and B one should choose any cubic roots from 
the corresponding expressions, which fulfill the following 
equation: 

AB = -p/3. (B5) 



In case this equation has real coefficients one should take 
the real values of the roots (if possible) . 
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